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Abstract.  This  paper  defines  a  program  realization  for  finding  the  optimal  geometry  of 
a  planar  Z-mechanism.  The  paper  shows  the  mathematical  procedure  for  defining  the 
objective  function,  the  constraint  function  and  the  search  fields  which  are  used  for 
solving  the  optimization  problem.  Starting  from  the  solutions  in  practice,  a  numerical 
example  is  given  to  determine  an  optimum  design  with  four  optimization  parameters. 
All  the  optimization  parameters  are  geometric  on  the  mechanism  for  defining  the  bucket 
position.  The  problem  is  solved  by  two  different  numerical  methods  -  the  method  of 
formal  search  of  hyperspace  (passive  scanning  method)  and  the  approximate  method  of 
Sequential  Quadratic  Programming  -  SQP  (applying  the  function  fminmax  of  the 
Matlab  optimization  toolbox).  The  verification  is  performed  with  animation  using 
software  for  geometric  modeling.  The  results  are  graphically  illustrated. 
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1.  Introduction 

For  working  with  bulk  material  the  wheel  loaders  have  a  mechanism  that  enables  the 
position  adjustment  -  an  inclination  angle  of  the  bucket.  At  lifting,  the  bucket  should  keep 
an  upward  position  without  changing  the  inclination.  This  can  be  achieved  by  adjusting 
the  geometry  of  the  parts  of  the  bucket  movement  mechanism.  One  of  the  mechanisms 
providing  for  it  is  known  as  the  Z  -  mechanism  (Fig.  1).  The  task  of  finding  a  solution  for 
the  special  mechanism  providing  for  a  special  request  is  known  as  the  optimal  synthesis 
according  to  the  Theory  of  Mechanisms. 

Solving  this  problem  is  a  common  example  of  research  and  development.  The 
similarity  of  different  plane  mechanisms  of  the  machines  offers  the  possibility  of  using 
the  same  optimization  procedure.  Nowadays  there  is  a  variety  of  mathematical  methods 
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in  use,  from  the  classical  one,  based  on  the  authorized  procedures  of  differential 
programming,  gradient  methods,  the  method  of  polyhedra  (Nelder-Mead  methods)  [1], 
approximate  methods  of  programming  and  sequential  quadratic  programming  [2],  up  to 
the  modern  methods  based  on  the  techniques  of  genetic  algorithm  [3],  ACA  algorithm 
(Ant  Colony  Algorithm)  [4]  and  simulation  [5], 

This  paper  first  presents  the  method  of  formal  searching  for  optimal  solutions  based 
on  the  computer  changes  of  the  geometric  parameters  of  the  mechanism,  followed  by  the 
method  of  Sequential  Quadratic  Programming.  Changing  the  parameters  in  the  method  of 
formal  searching  is  performed  incrementally  by  taking  discrete  values  of  independent 
parameters  in  an  increasing  order  (scanning  technique)  [6],  [7]. For  the  evaluation  of  the 
effectiveness  of  the  passive  formal  search  method,  the  method  of  Sequential  Quadratic 
Programming  (SQP)  is  applied  in  this  paper. 

The  complexity  of  the  general  approach  to  the  problem  can  be  seen  in  the  work  by 
Yang  Yu  [8]  where  the  optimization  is  performed  with  twenty-two  parameters.  Besides 
geometric  optimization  parameters,  the  synthesis  can  also  introduce  the  constraints  of 
internal  mechanism  forces.  Such  an  approach  can  be  found  in  the  work  of  Gao  Xiuhua  [9] 
(derived  using  the  Adams  software  package),  and  Rongyi  Zhang  [10],  using  the  Matlab 
programming  environment  with  the  use  of  a  combination  of  three  optimization  algorithms. 
Because  of  the  complexity  of  the  general  approach,  this  paper  deals  with  the  mechanism 
synthesis  based  on  a  small  number  of  optimization  parameters. 


2.  Geometric  Model  of  the  Z-Mechanism 


A  planar  mechanism  of  a  manipulator  loader  consists  of  an  arm  (1)  with  a  bucket  (4) 
and  a  lifting  hydraulic  actuator  (5),  Fig.  1.  The  arm  (1)  possesses  a  Mechanism  for 
Defining  the  Bucket  Position  (MDBP)  that  does  loading,  transport  position  holding  and 
emptying  of  the  bucket.  This  mechanism  consists  of  a  tilting  hydraulic  actuator  (6),  a 
two-arm  lever  (2)  and  a  coupling  rod  (3)  (Fig.  1). 

Solving  the  optimization  problem  involves 
finding  a  geometry  of  mechanism  for  defining  the 
bucket  position  to  ensure  minimal  change  in  the  tilt 
of  the  bucket  when  the  arm  of  the  loader  is  lifting 
from  the  minimum  to  the  maximum  slope  angle  /?, 
min  {FC}.  From  there,  it  is  possible  to  search  for 
optimization  parameters  in  the  MDBP  mechanism 
geometry.  The  mechanism  for  the  lifting  arm  is 
designed  from  the  given  operating  conditions  of  the 
bucket  position  and  is  not  subject  to  optimization. 
The  optimization  parameters  can  therefore  be  the 
position  of  the  bearing  of  the  hydraulic  actuator  (6), 
the  geometry  of  the  two-arm  lever  (2)  and  the 
length  of  the  coupling  rod  (3).  In  this  study,  four 
parameters  are  chosen  to  optimize  the  MDBP 
mechanism  related  to  the  position  of  the  bearing  of  the  hydraulic  actuator  (6)  and  linkage 
O23  on  the  lever  (2). 


Fig.  1  Geometric  model  of  the 

Z  -  mechanism  for  defining 
the  bucket  position  (MDBP) 
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For  an  efficient  solution  of  multidimensional  problems  with  constraints,  numerical 
methods  with  passive  searching  can  be  used  [6].  These  methods  use  a  fixed  interval  of 
indetermination  (H),  where  the  field  of  change  of  independent  variables  z  =  {zl,z2,z3,z4} 
is  divided  by  the  selected  number  of  points  equal  to  the  number  of  sub -areas.  These  points 
serve  as  the  basis  for  the  construction  of  n  dimensional  space.  The  objective  function 
FC(z)  is  calculated  at  each  node  of  the  grid  thus  formed.  With  the  analysis  of  all  solutions 
the  one  that  is  the  most  advantageous  is  chosen.  The  advantage  of  the  passive  search 
method  is  in  the  reliability  of  searching.  A  disadvantage  of  multidimensional  problems  of 
passive  searching  lies  in  the  number  of  calculation  procedures.  Therefore,  this  paper 
employs  the  procedure  of  choice  (constriction)  of  the  field  of  search  based  on  the  testing  of 
the  objective  function.  The  testing  determines  the  solution  with  the  good  initial  values  of 
FC(z)  based  on  a  large  interval  of  indetermination.  After  identifying  the  "good  areas"  of 
solution,  the  subject  of  solution  is  the  grid  with  a  fine  interval  of  indetermination.  Such  a 
model  with  a  reduced  number  of  optimization  parameters  (zi,  z2,  z3,  z4)  and  numerical 
implementation  can  give  sub-optimal  solutions  thanks  to  the  simplicity  of  the  algorithm 
search.  The  synthesis  procedure  for  an  MDBP  mechanism  is  illustrated  in  Fig.  1. 


3.  Mathematical  Formulation  of  the  Problem 

The  formation  of  a  mathematical  model  for  optimization  according  to  Fig.  1  includes 
defining  a  planar  geometric  model  of  the  manipulator.  The  manipulator  is  a  planar  (2D) 
mechanism  with  an  arm  for  carrying  the  bucket  and  a  mechanism  for  defining  the  bucket 
position.  A  mechanism  for  defining  the  bucket  position  is  labeled  with  darker  shades  in 
Fig.  1.  The  whole  mathematical  model  uses  the  global  Cartesian  x-y  coordinate  system 
and  local  coordinate  systems  of  the  parts  of  the  mechanism  shown  in  Fig.  2.  Fig.  2  shows 
four  parts  of  the  mechanism.  The  planar  form  of  geometry  of  these  two  coupled 
mechanisms  with  the  selected  marked  positions  of  parts  is  given  in  Fig.  3. 


Fig.  2  Parts  of  MDBP:  a)  the  arm  of  the  manipulator,  b)  the  two -arm  lever, 
c)  the  coupling  rod,  d)  the  bucket 
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Let  us  define  the  geometry  of  the  rigid  mechanical  system  in  Fig.  l:The  coordinates 
of  linkage  On  are  given  values  Xn,  Yu  of  some  type  of  a  loader.  The  coordinates  of 
linkage  0|2  at  the  top  of  the  arm  are  computed  as: 

xn  =xn  +a2 -cos/?  +  a3 -cos(90-/?)  (1) 
Yl2  =Yu+a2- sin /3-a3cos/3  (2) 

The  coordinates  of  bearings  of  the  two-arm  lever,  the  position  of  linkage  022  are 
computed  as: 

X22  =  Xu       -cos/?  (3) 

Y22  =YU  +a,  -sin ft  (4) 

The  coordinates  of  linkage  (92i  are  determined  by  an  intersection  of  the  two  circles 
with  the  center  in  linkage  Ob  whose  radius  is  c6  -  the  length  of  the  hydraulic  actuator  (6) 
and  the  circle  with  the  center  in  linkage  <922  whose  radius  is  equal  to  the  length  of  the  arm 
/,  of  the  lever  (2).  We  consider  that  the  length  of  the  hydraulic  actuator  (6)  is  constant. 
The  coordinates  of  linkage  <92i  are  obtained  by  solving  the  equation: 


Fig.  3  Geometric  model  of  the  manipulator  of  the  loader 


The  distance  between  the  centers  of  the  circles  is  computed  as: 

D  =  ^(Xb-X22f+{Yb-Y22f 


(6) 
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Value  F  and  the  coordinates  of  linkage  02i  are  computed  as: 

F  =  ^(D  +  ll+c6)-(D  +  l1-c6)-(D-l,+c6)-(ll-D  +  c6)  (7) 
r       =  X22  +X6  +(X6  -X22)-(l^ -cl)  ±2Y^F 

21(1,2)  2  2-Z)2  ^2 


^21(1,2,=— ^  +  ^3  +2         n2  W 


Let  us  determine  the  coordinates  of  linkage  (923-  Angle  ^ ;  between  direction  r  and  d  of 
the  two-arm  lever  is  computed  as: 

L 

rii=arctg——  (10) 
l2  /j 

Variable  angle  tj  (between  the  arm  and  the  two-arm  lever),  if  is  X22  -  X2\  >  0  is  then 
computed  as: 

n  =  arctg^~_Y*  (11) 

A22  A21 

If  X22  -  X2\  =  0,  angle  rj  =  90°,  and  if  X22  -  Xa<0,  angle  tj  is  computed  as: 

7  =  180° -arctg  J21  (12) 

21  22 

Length  d  of  the  arm  of  the  two-arm  lever  (Fig.  2b)  is  computed  as: 

d  =  J(l2-llf  +  %  (13) 
Coordinates  of  linkage  O22,  are  computed  as: 

X23  =  X22  +d  cos(77  +  771)  (14) 

Y23  =  F22 -<i  sin(77  +77j)  (15) 

The  introduced  optimization  parameters  zh  z2,  z3,  z4  correspond  to  the  required 
geometry  of  the  Z  -  mechanism  by  the  following  relations: 

zx=X6  (16) 

z2=Y6  (17) 

z3=/2-A  (18) 

z4=h  (19) 
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The  coordinates  of  linkage  04  (X4,  Y4)  are  determined  using  the  intersection  of  the 
circle  with  the  center  of  linkage  (923  whose,.  The  solution  of  coordinates  X4,  Y4  is  obtained 
from  the  equation:  radius  is  the  length  of  coupling  rod  C\  and  the  circle  with  the  center  of 
linkage  On  whose  radius  is  length  b: 


(X4-X23)2+(Y4-Y23f=cf 
(X4-Xi2)2+(Y4-Yl2)2  =b2 


(20) 


The  distance  between  the  centers  of  the  circles  and  value  Fj  are  computed  as: 

Di  =  ^Xn-X23)2+(Y23-Yu)2  (21) 


F^-^P,  +c,  +b1)-(Dl  +c,  -^)-(A  -q  +bl)-(c1  -D  +  b,) 


(22) 


Coordinates  of  linkage  04 : 


Xl2+X23    (Xl2-X23)-(cf-b2)  Y23-Yl2 
^=—2  —  +  ^  ±2— -Fl  (23) 

= ^+ {Yu  ~Yf£  ""'2)  +2^^  (24) 

Bucket  tilt  angle  for  the  current  value  of  arm  slope    is  computed  as: 

a{P)  =  arctg  J4~^2  (25) 

The  back  angle  of  bucket  y  is  determined  by  the  position  of  the  angle  of  bucket  a  and 
constant  a4  (Fig.  2d): 

y  =  a  -  a4  (26) 

The  objective  function  is  the  minimum  change  in  the  slope  of  the  bucket  while  the 
arm  of  the  loader  is  lifting: 

FC{zl  ,z2,z3,z4)  =  amm-  a„  (27) 

The  constraint  functions  are  determined  by  inequalities: 

Gl=zl-A>0  G5=z3-E>0 
G2=z,-B<0      G6  =  z3-F<0 

2         1  6         3  (2g) 

G3  =  z2-C>0  G7=z4-U>0 
G4=z2-D<0  Gs=z4-H<0 

where:  A,  B,  C,  D,  E,  F,  U,H-  are  the  constants  of  the  search  space. 
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4.  Algorithm  and  Program 

Fig.  4  presents  the  algorithm  of  global  numerical  procedures  developed  for  the  solution 
of  the  optimal  mechanism  synthesis.  The  algorithm  shows  that  the  scanning  method  is 
realized  by  an  increment  change  of  the  optimization  parameters  in  a  permissible  four- 
dimensional  hyperspace.  The  permissible  hyperspace  is  empirically  determined.  The  initial 
setting  of  geometry  of  both  mechanisms  is  provided  through  the  files  DAT1-DAT6.  The 
solutions  of  objective  function  FC  for  discrete  values  of  independent  parameters  z\  -f  Za  are 
entered  in  the  file  DAT7.  By  comparing  these  solutions,  the  solution  that  is  characterized  by 
the  minimum  value  of  objective  function  FC  is  chosen.  Then  the  obtained  solution  is 
checked  from  the  point  of  numerical  accuracy,  satisfying  the  given  nominal  characteristics 
of  the  reach  and  lift  height.  After  that,  the  continuity  of  coupling  of  the  parts  of  mechanism 
for  defining  the  bucket  position  is  also  checked.  The  following  numerical  example 
illustrates  the  application  of  this  mathematical  model,  the  efficiency  of  the  algorithm  and 
the  practical  realization  of  optimal  synthesis. 


C 


START 

ZJZ. 

Programme  title 


j 


-_PATL 

~datT 


DAT3  . 


1  1 

Introduction  of 

initial  data  about 

^~LT  DAT5 
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Setting  parameter 
optimization 
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Fig.  4  The  algorithm  of  the  program  for  optimal  numerical  synthesis  of  MDBP 
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5.  Numerical  Example 

The  aim  of  this  study  is  to  search  for  the  design  of  the  mechanism  for  defining  the 
bucket  position  which  minimizes  the  angle  of  the  bucket  tilt  during  the  bucket  lifting 
operation.  The  initial  solution  is  based  on  freely  assumed  geometry  which  is  in  this  case 
taken  from  photographs  and  sketches  [11]  (Komatsu  WA320).  Four  parameters  are 
selected  to  optimize  the  synthesis:  x,  y  coordinates  of  bearings  of  the  hydraulic 
mechanism  for  defining  the  bucket  position  and  x,  y  local  coordinates  of  linkage  023  of 
the  two-arm  lever.  The  procedure  for  the  algorithm  program  flow  is  fast  and  numerically 
stable.  The  stability  is  based  on  the  simplicity  of  analytical  modeling  geometry. 

Angle  P  is  the  slope  angle  of  the  arm  in  relation  to  the  horizontals  and  ranging,  in  this  case, 
from  -24°  to  62°.  Parameters  zb  z2  represent  the  search  field  of  bearings  position  (96,  where  the 
geometry  is  in  the  following  ranges:  in  the  x-axis  direction  Zimin=0-94  m;  Zimax=1.14  m 
(Az!=0,2  m)  (A=0.94  m;  B=1.14  m)  and  in  the  y-axis  direction  z2min=1.645  m;  z2max=1.845  m 
(Az2=0.2  m)  (C=1.645  m;  D=  1.845  m).  The  allowable  field  defined  by  constants  A,  B,  C,  D  is 
chosen  on  the  basis  of  the  available  space  on  the  arm.  By  choosing  M  =  50  points  of  change  of 
one  independent  parameter  of  optimization,  the  interval  of  indetermination  is: 

R  —  A       1  1 4  —  D  Q4 
H,=2— -  =  2  '        '     =0,00816  [m]  (29) 
1      M-l  50-1 

Parameters  z3,  z4  represent  the  search  field  of  the  position  of  linkage  <923  in  the  local 
coordinate  system  of  the  two-arm  lever,  the  size  of  change  is  in  the  area:  z3mjn=1.4  m; 
Z3max=l-5  m  (Az3=Q.\  m),  (£=1.4  m;  F=1.5  m)  in  the  x-axis  direction,  and  z4mjn=0.175  m; 
Z4max=0-275  m  (Az4=Q.\  m),  (J/=0.175  m;  H=Q.215  m)  in  the  y-axis  direction.  The  field  is 
empirically  chosen. 

Initial  zeros  are:  z,=1.04  m,  z2=1.745  m,  z3=1.45  m,  z4=0.225  m.  These  areas  of 
solutions  of  the  optimization  task  determined  the  allowable  space  of  search-scanning. 
The  dividing  of  the  permissible  field  of  parameters  z3  and  z4  by  the  N  =  50  points  is 
selected,  as  well  as  the  dividing  of  the  permissible  field  of  parameters  z\  and  z2  by  M  =  50 
points.  In  this  way  a  number  of  combinations  of  K=M  2,N  2=6250000  are  obtained.  From 
the  file  DAT-8  we  can  read  the  optimal  parameters  of  the  mechanism  for  defining  the 
bucket  position.  Based  on  these  parameters,  it  is  possible  to  draw  curves  of  objective 
function  FC,  depending  on  the  angle  of  the  slope  of  the  arm  /?  for  the  optimal  and  initial 
solution  (Fig.  5).  Starting  from  the  fixed  size  of  the  loader  mechanism  given  in  Table  1, 
the  MDBP  optimal  solutions  are  obtained,  given  in  Table  2.  The  curves  of  the  objective 
function  of  the  required  MDBP  are  plotted  for  the  optimal  solution. 


Table  1  Initial  data  of  model 


a  i  (m)  a2(m) 

a3(m)        h(m)         l2(m)  Z3(m) 

fli  (m)         b\  (m) 

1.800  2.550 

0.865         0.720         1.450  0.225 

0.765  0.370 

Table  2  Values  of  the  obtained  optimal  solution 

Optimal  parameters 

Objective  function 

zi  (m)        z2(m)      z3(m)  z4(m) 

FCC) 

Initial  zeros 
Optimal  values 

1.04          1.745      1.45  0.225 
0.973        1.678      1.400  0.175 

7.364 
1.707 
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Fig.  5  Curve  of  initial  FCo  and  curve  of  optimal  solution  FQ 

The  verification  is  performed  by  simulation  of  the  working  mode  (meshing)  of  the 
mechanism,  in  the  SolidWorks  software,  by  using  the  option  blocks  which  is  suitable  for 
the  simulation  of  planar  mechanisms.  This  achieved  the  verification  of  the  behavior  of  the 
mechanism  in  successive  positions,  Fig.  6.  It  provides  a  functional  control  study  of  the 
mechanism  and  eliminates  errors. 


L 


Fig.  6  Simulation  of  the  obtained  optimal  solution 


6.  Minimization  Using  SQP  Method 

In  order  to  assess  the  benefits  from  the  chosen  method  of  minimizing  the  objective 
function,  the  method  of  the  Sequential  Quadratic  Programming  (SQP)  is  used.  Minimization 
is  performed  by  programming  the  objective  function  and  constraint  function  in  Matlab  by 
using  the  program  function  (command)  fminimax  as  part  of  the  Madab  optimization  toolbox. 
The  general  quadratic  programming  problem  is  defined  as: 
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min  f(z)=cTz+-zTH-z  (30) 

Where:  A  ■  z  <  b  ,  z  >  0  (31)  are  the  constraints  in  the  form  of  inequality, 
c  -    the  n-dimensional  vector  that  describes  the  coefficients  of  linear  terms  of  the 

objective  function, 
z  -    the  n-dimensional  vector  of  unknown  values, 

H—  the  (nxn)  symmetric  matrix  that  describes  the  coefficients  of  square  terms  of  the 

objective  function, 
A  -  the  (to  x  n)  matrix  of  constraints,  and, 
b  -   the  n-dimensional  vector  of  constraints  [12]. 

5.1.  Function  minimax 

Function  minimax  (fminimax)  is  part  of  the  Optimization  Toolbox  in  Matlab,  which 
finds  the  minimum  of  problems  determined  by  the  function: 

c(z)  <  0 
ceq(z)  =  0 
A-z<b  (32) 
Aeq ■ z  = beq 
lb<z<ub 

Where  the  variables  on  the  right  side  are  the  constraints: 
c(z),  ceq(z)  -  vectors  of  constraints  in  the  form  of  inequality  and  equality, 
A  -  matrix  and  b  -  vector  are,  respectively,  the  coefficients  of  constraints  in  the  form  of 
linear  inequalities, 

Aeq  -  matrix  and  beq  -  vector  are,  respectively,  the  coefficients  of  linear  constraints  in  the 

form  of  equity, 
lb  -  vector  or  matrix  of  lower  bounds  of  constraints, 
ub  -  vector  or  matrix  of  upper  bounds  of  constraints. 

In  the  model  that  is  described  by  equations  1-28,  we  use  only  constraints  lb,  ub  -  the 
vectors  of  lower  and  upper  bounds  of  the  search  field  of  the  optimization  parameters 
solution  (A,  B,  C,  D,  E,  F,  U,  H),  whose  values  are  listed  in  the  previous  numerical 
example  [13]. 

Based  on  the  function  fminimax,  where  it  is  necessary  to  define  the  objective  function 
and  the  constraints,  the  desired  optimal  parameters  are  obtained: 

z  =  fminimax  ( fun,  z0,  A,  b,  Aeq,  beq,  lb,  ub)  (33) 

where:  z  -  required  parameters,  fun  -  objective  function,  zo  -  initial  parameters  A,  b,  Aeq, 
beq,  lb,  ub  -  constraints.  These  values  are  defined  above  in  equations  1-h28,  since  in  this 
case  constraints  A,  b,  Aeq  and  beq  are  not  used.  They  are  written  as  an  "empty"  matrix. 
The  inputs  of  constraints  are  in  the  form  of  the  lower  bounds  -  lb  and  the  upper  bounds  - 
ub.  In  this  case,  the  line  of  the  code  of  function  fminimax  looks  like: 

parameter s= fminimax ( ' Ob j  ectf un ' , parametriO ,  [],[],[],[], lb , ub)  (34) 


min  max  Ft  (z)  such  that  is  < 

z  i 
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SQP  has  performed  the  verification  of  the  previous  procedure  and  the  optimization 
achieved  by  using  the  passive  formal  search.  The  results  obtained  by  the  SQP  method  are 
shown  in  the  following  diagram,  Fig.  7,  and  a  table  of  optimal  solutions. 

Table  3  Values  of  obtained  optimal  solutions  by  using  the  function  fminimax 


Optimal  parameters 

Objective  function 

Zi  (m) 

z2(m)  z3(m) 

Z4(m) 

FCC) 

Initial  zeros 

1.04 

1.745  1.45 

0.225 

7.364 

Optimal  solution 

1.0176 

1.7188  1.4095 

0.2002 

0.3518 

SOP 

method 
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Fig.  7  Curve  of  initial  FC0  and  optimal  FQsolutions,  obtained  by  using  the  Minimax  function 


7.  Conclusion 

1.  The  mathematical  model  of  optimization  operates  with  four  independent 
parameters  out  of  the  potential  eight.  Even  four  parameters  of  optimization  give  a 
solution  of  objective  function  (change  of  the  tilt  is  less  than  two  degrees)  good 
enough  for  practical  performance.  By  increasing  the  number  of  independent 
parameters  to  possible  eight,  the  quality  of  the  optimal  solution  can  be  increased, 
with  a  complex  mathematical  procedure. 

2.  It  is  recommended  to  use  an  optimal  model  of  the  MDBP  mechanism  for  factory 
design  because  it  is  an  important  design  requirement. 

3.  Improving  the  model  can  be  done  by  introducing  the  constraint  function  of 
internal  forces  in  the  MDBP  which  would  require  the  extension  of  this  task  in  the 
domain  of  static  analysis. 

4.  The  SQP  procedure  has  given  a  slightly  more  accurate  solution  which  is  a  result  of 
a  freely  adjustable  step.  The  formal  procedure  limits  the  accuracy  (quality  of 
solutions)  by  selecting  the  interval  of  indetermination.  The  quality  of  a  solution 
can  be  improved  by  increasing  the  number  of  search  combinations. 

5.  The  computing  time  of  implementation  in  each  case  is  less  than  one  second. 
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OPTIMALNA  SINTEZA  MANIPULATORA 
KORISTECI  DVE  UPOREDNE  METODE 

U  radu  je  data  programska  realizacija  za  trazenje  optimalne  geometrije  ravanskog  Z- 
mehanizma.  Rad  pokazuje  matematicku  proceduru  definisanja  funkcije  cilja,  funkcija  ogranicenja, 
oblasti  pretrazivanja  kojima  se  resava  zadatak  optimizacije.  Polazeci  od  resenja  u  praksi,  dat  je 
numericki  primer  odredivanja  optimalnog  dizajna  sa  cetiri  parametra  optimizacije.  Svi  parametri 
optimizacije  su  geometrijski  na  mehanizmu  za  odredivanje  nagiba  kasike.  Zadatak  je  resavan 
dvema  razlicitim,  numericki  metodama  -  metodom  formalnog  pretrazivanja  hiperprostora  (metoda 
pasivnog  skeniranja)  i  aproksimativnom  metodom  kvadratnog  sekvencijalnog  programiranja  -  SQP 
(primenom  fininmax  funkcije  iz  Matlabovog  optimizacionog  toolbox-a).  Verifikacija  resenja  je 
vrsena  animacijom  u  programu  za  geometrijsko  modeliranje.  Rezultati  su  graficki  ilustrovani. 

Kljucne  reci:  Pasivno  trazenje,  Z  -  mehanizam,  utovarivac,  rudarska  masina,  optimizacija 
mehanizma,  numericko  resavanje,  SQP 


